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I. INTRODUCTION 

This paper focuses on the quantum gravity extension of the inflationary scenario provided 
by loop quantum cosmology (LQC) [1-4]. We follow the framework introduced by Agullo, 
Ashtekar and Nelson in a series of papers [5-7]. (For other approaches to the early universe 
in loop quantum cosmology see [8-24].) In short, in this framework quantum gravitational 
effects dominate the Planck era of the universe, and a quantum bounce appears replacing 
the classical big bang singularity. For dehniteness, the matter content of the universe is 
assumed to consist of a massive scalar (“inflaton”) held, although other forms of the potential 
U(0) can be accommodated without altering the conclusions. Shortly after the bounce, 
quantum gravity effects gradually lessen and the potential energy U(0) begins to prevail. 
This potential-dominated phase brings the universe, under quite generic circumstances, to a 
phase of inhation [25, 26]. Therefore, LQC provides an interesting arena to incorporate the 
highest energy density and curvature stages of the universe into cosmological models, where 
questions about Planck-scale physics and initial conditions for inhation can be addressed 
squarely. 

The theory of cosmological perturbations in the Planck era was developed in [6], fol¬ 
lowing previous results in [27]. It was then applied in [7] to show that, generically, the 
pre-inhat ionary evolution makes scalar and tensor perturbations reach the onset of inhation 
in an excited state, rather than in the Bunch-Davies vacuum often assumed in the cosmol¬ 
ogy literature. Consequently, the inhationary primordial spectra that source the cosmic 
microwave background (CMB) anisotropies acquire some extra features with quantum grav- 
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itational origin. As described in detail in [7] and summarized in section II, the effects of the 
pre-inflat ionary evolution are more important for infrared modes (which correspond to large 
angles in the CMB), and the size of these effects depends on the parameters of the model; 
specihcally the inflaton mass m and the “initial” value at some reference time—for which 
we choose the bounce time—of the inflaton field, 0(tBounce) = 0 b- 

The hrst phenomenological exploration of the spectrum of scalar and tensor perturbations 
under this framework appeared in [7]. For phenomenological interest, this analysis focused 
on a small region of the parameter space of the theory, consisting of values of 0b near to 
its minimum possible value and the inflaton mass m that is commonly used in inflationary 
cosmology. The main goal of this paper is to extend the study to the full parameter space 
spanning the plane (0b, m). 

The goals are multiple: (i) compute the power spectra of scalar and tensor perturba¬ 
tions for the theoretically allowed values of 0b and m; ( 2 ) identify which region of the 
(0B,'m) plane passes all observational constraints; (<?) localize the subspace which, in ad¬ 
dition to being observationally allowed, contains significant corrections, originating in the 
pre-inflat ionary evolution of LQC, to the standard inflationary picture; and {4) analyze in 
detail the predicted observational signatures, which will mostly involve tensor modes. 

The main challenge of the analysis presented here is computational. Time—and 
memory—intensive computations using high-performance computing are required to explore 
the most interesting region this space. But this effort is necessary to understand completely 
the predictions of the theory across the parameter space. 

Our results are in agreement with those obtained previously in [7] when we restrict to 
the range of parameters explored there, but interesting new hndings arise in other regions 
of the parameter space. We summarize here the most important points: 

• Two scales dictate the form of the LQC-corrected power spectrum; First, /clqc, the 
momentum scale associated with the spacetime curvature at the bounce, where it 
attains its maximum value; and second, ki, the momentum scale associated with the 
spacetime curvature at the onset of inflation. (It always happens that /clqc 3> ki.) 

• The LQC power spectrum is oscillatory and its average has an amplitude that is 
amplihed with respect to the standard predictions of slow-roll inflation for modes 
h ^ k < /clqc, but is in agreement with them otherwise (see Fig. 2). 

• The present values of the physical scales /ci/oq and /clqc/uo (uq is the scale factor 
today) depend on the amount of expansion that has occurred since the time of the 
bounce, which in turn depends on the values of the parameters 0b and m. Therefore 
0B and m control whether the modes affected by the pre-inflationary dynamics of LQC 
fall within the window of modes which are observable today. 

• The region of parameter space which is observationally viable is 

0.8 < 0B 

1.1 X 10“® < m < 1.5 X 10“®. 


( 1 . 1 ) 


• The region which, in addition to being observationally allowed, also shows non- 
negligible LQC modihcation is an approximately one-dimensional subset of (1.1), 
which satisfies 0 b ~ 1-3 x 10“®/m (see Fig. 4). 
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• LQC corrections to the power spectra tend to make the tensor spectral index rit more 
negative and to produce a positive running ag of the scalar spectral index. 

• The corrections tend to reduce the the tensor-to-scalar ratio r, which serves to alleviate 
the observational constraints on the potential. 

• The LQC corrections modify the inflationary consistency relation r/rit = —8. 

• The particular choice of initial data for quantum scalar and tensor perturbations has 
very little impact on the above conclusions, at least for the reasonable choices of initial 
vacuum state that we have considered in this paper. 

In the rest of this paper we provide the details and summarize the calculations leading 
to these conclusions. We begin in section II by summarizing the pre-inflationary physics 
in LQC, both for the background spacetime and for cosmological perturbations, and by 
describing the main features of the resulting power spectra for a typical evolution. In 
section III we report the results of exploring the predictions of LQC across the parameter 
space and their relation with observations. Section IV analyzes the sensitivity of the results 
to the initial quantum state of scalar and tensor perturbations. In section V we discuss our 
results and add some hnal comments. 

Throughout this paper all numerical values are given in Planck units, in which c = 
G = h = 1. Consequently the Planck length, time, and mass all equal unity: £pi = tp\ = 
mp\ (= \Jhc/G) = 1. However, we will retain G and h explicitly in theoretical expressions 
to emphasize their physical content. 

Remark: In previous analysis [5, 7] WMAP’s 7-year observational data was used to 
compare with observations and to £x some parameters—for instance the inflaton mass. In 
this paper we use the more recent 2013 Planck results [28] {Planck 2015 results [29] are very 
similar). Therefore, some care is needed in directly comparing numerical values across the 
two analyses. 


II. SUMMARY OF THE PRE-INFLATIONARY EVOLUTION OF LOOP 

QUANTUM COSMOLOGY 

The goal of this section is twofold. First, we summarize the physics of the spacetime 
evolution of the early universe in loop quantum cosmology, as well as the equations gov¬ 
erning the evolution of first-order scalar and tensor perturbations thereon. The material 
has been explained in full detail in the original references [6, 7, 25, 26] and summarized in 
review articles [2, 25, 30]; we therefore provide here only a succinct summary and refer the 
reader to those references for further subtleties and details. Second, we analyze the main 
features of the power spectra for a representative concrete choice of initial parameters. This 
analysis will provide an understanding of the origin and characteristics of the corrections 
that LQC introduces to the primordial spectra. While some of this material was computed 
and analyzed in detail in [7] as well, we also discuss a number of new features of interest. 
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A. Dynamics of the FLRW spacetimes in LQC 


Since this paper focuses on results directly relevant for observations, we restrict to spa¬ 
tially flat, FLRW spacetimes. The gravitational held is sourced by a single scalar held 0, the 
inhaton, with an ehective potential V (0) that we choose to have the simple quadratic form 
R(0) = As mentioned previously, other choices are certainly possible and, although 

the concrete numerical values obtained below would change, our hndings would remain 
qualitatively unaltered; the underlying reason for this is that the LQC ehects on primordial 
perturbations originate from quantum-gravitational ehects that are largely insensitive to the 
specihc form of V{(j)). 

In LQC the quantum homogeneous and isotropic spacetime is described by a quantum 
state To (a, 0) which is a complex function of the classical scale factor a (or equivalently its 
cube, V = a^, which is more commonly used in the LQC literature) and the homogeneous 
part of the inhaton held 0. Among the many states To(a, 0) contained in the LQC Hilbert 
space [31-34], of particular physical relevance are those which are sharply peaked around 
a classical trajectory at late times, when the curvature of the universe is well below the 
Planck scale. Most of these states continue to be sharply peaked during the entire dynamical 
trajectory [32, 33], including the Planck era. More importantly, the evolution of the peak 
of the wavefunction of such states can be accurately described by ehective equations that 
very much resemble the equations of general relativity, apart from some extra terms that 
account for the quantum gravity corrections [see Eqs. (2.1) and (2.2) below]. One can 
therefore explore the phenomenological consequences of the theory without the need to 
solve the complicated discrete equations of the full theory. For this reason, the ehective 
dynamics have been used in essentially all phenomenological analysis performed so far in 
the literature. We do the same in this paper. ^ 

The ehective equations of LQC were derived in [36-39]. The modihed Friedmann equa¬ 
tions is given by 

\ Pmax / 

the modihed Raychaudhuri equation is 


( 2 . 1 ) 


a 

a 


AttG 


p 1-4- 


P 


Pn 


-AttGPI 1 



( 2 . 2 ) 


and the equation of motion for the inhaton held, which takes the same form as in the classical 
theory, is 

dV 

+ 3H + — = 0 . (2.3) 

d0 

In these equations p = |0^ -|- V (0) and P = |0^ — V (0) are the energy and pressure density, 
respectively, of the scalar held; p^ax = 0.41ppi is the upper bound of the energy density in 


^ One can, however, legitimately argue that other quantum states 5'o(a, 0) showing large quantum fluctu¬ 
ations in the Planck era may be of physical interest. States ikolttid) which are not sharply peaked are 
not accurately described by the effective equations, and to explore their phenomenology one is forced to 
face the full discrete equations of LQC, which are significantly more complicated. The analysis that we 
present in this paper has been extended to such more generic family of quantum states, and the results 
will be presented in a separate publication [35]. 
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LQC, where ppi is the Planck energy density. Note that pmax is proportional to thus in 
the classical limit p^ax diverges and Eqs. (2.1) and (2.2) reduce to the classical Friedmann 
and Raychaudhuri equations. 

The relevant parameter spaee. The previous equations can be solved numerically by 
specifying initial data for a, a, 0 and 0 at a given time, together with a value of the inflaton 
mass m. The bounce is a convenient time to specify initial data because the universality of 
the solutions there reduces the number of free parameters, as we now explain. First, note 
that in a flat FLRW spacetime only ratios between values of the scalar factor at different 
times have objective physical meaning, not the value of a{t) itself at any one time. We are 
therefore free to rescale a at our convenience. A convenient choice is a{t-B) = 1, at the bounce 
time t-Q- Second, note that, because the potential R(0) is symmetric, the transformation 
(0(fB), 0(tB)) —> (—0(tB), —0(tB)) does uot alter the physics. We can therefore restrict 
our solutions to 0(tB) > 0 without loss of generality. Third, at the bounce we always have 
a = 0. Fourth, from Fq. (2.1) it is easy to see that the value of p at the bounce must equal 
Pmax- This implies that 0(tB) and 0(tB) are related by pit-o) = |0(tB)^ + ^?7i^0(tB)^ = Pmax- 

In summary, solutions to Fqs. (2.1), (2.2) and (2.3) form a two-parameter family labeled 
by couples (0(fB), "nr). (From now on we will denote 0(tB) by 0 b.) Because the energy density 
is bounded above, the parameters |0b| and m are as well—they must satisfy m|0B| < 0.90. 
For dehniteness, we will assume 0 b > 0 because the sign will not make any qualitative 
difference in our analysis. We conclude, therefore, that the relevant parameter space for this 
paper is the set of couples (0b, nr) satisfying 


0 < m0B < 0.90 . 


(2.4) 


We now briefly summarize the main qualitative features of the solutions to Fqs. (2.1), 
(2.2) and (2.3), which are illustrated in Fig. 1. All solutions experience a bounce at which 
H = t) and p = Pmax- The bounce has a quantum gravitational origin that makes it 
independent of the matter content of the theory, and in particular independent of the form 
of the potential V (0). The effective spacetime geometry is symmetric around ts for a kinetic- 
dominated bounce (0(tB)^ ^ m^0|), which, as explained in section III, turns out to be the 
most interesting regime for possible new predictions. We will therefore focus our discussion 
on this regime. Immediately after the bounce, the Hubble rate JT grows from zero to its 
maximum value TTmax ~ 0.93, which is attained around 0.2 Planck seconds after the bounce. 
Because JT > 0, this period is commonly known as super-inflation (there is an essentially 
symmetric period of super-deflation prior to the bounce). After super-inflation, IT becomes 
negative as the inflaton held keeps climbing up the potential at the expense of its kinetic 
energy, gradually entering a potential-dominated regime. Around 10® Planck seconds after 
the bounce, the inhaton loses its last remaining kinetic energy, stops moving upward and 
begins a phase of slow-roll back down the potential; this is the onset of inhation. At this 
time the energy density has decreased approximately eleven orders of magnitude since the 
bounce, and consequently the quantum effects of gravity are negligible. The duration of the 
slow-roll phase depends on how high the inhaton has been able to climb up the potential; 
this grows monotonically with 0 b. 

In order to dehne more precisely the times at which the phase of slow-roll inhation begins 
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FIG. 1: A typical evolution of the energy density and Hubble rate of the background spacetime 
(in Planck units). Note that both axes are logarithmically scaled. The LQC modification to the 
classical Friedmann equation is evident near the bounce (since classically p = but by 2 

Planck seconds the behavior coincides with the classical trajectory. 

and ends, we introduce the first-order slow-roll parameters^ 

= -iL x = JL_ 

^ ’ 2HH 

Slow-roll inflation is said to have begun when both of these parameters are much smaller 
than unity; for concreteness, we define the onset of slow-roll inflation in this work as the 
time ton such that e(fon) = 0.1 and 6{ton) = 0.1, and the end of slow-roll inflation tend as the 
first time after ton when e(tend) = 1- 


(2.5) 


B. Evolution of cosmological perturbations 

In the standard inflationary scenario, cosmic non-uniformities are described by hrst-order 
scalar and tensor perturbations. Scalar perturbations can be conveniently described by the 
Mukhanov-Sasaki variable Q{x),^ and we will collectively denote by T{x) the two degrees 


^ The slow-roll parameters defined here, which are calculated in terms of the Hubble parameter, differ from 
another frequently-encountered set of slow-roll parameters, in terms of the potential V: 


C — 1 /dV\2 

tv — 2 vH 90 / > 


A — 1 9^V 
"V — V 902 • 


While these two sets of parameters are distinct, they are related if H and H are sufficiently small by 
tv ~ e, (5v ~ e — (5. 

^ When the matter sector is a single scalar field, this variable relates to the standard comoving curvature 
perturbation 7^ by 72. = ^Q. As explained in section V of [7], to analyze the pre-inflatioanry universe it 
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of freedom in tensor perturbations. Because the background energy density is well below 
the Planck scale during inflation, perturbations are accurately described as quantum helds 
propagating in the classical FLRW inflationary spacetime. This mathematical framework 
of quantum field theory in classical but curved spacetimes was developed in the late 1960s 
and ’70s, and since then it has been successfully applied to multiple interesting physical 
situations. But prior to inflation, and particularly in the vicinity of the bounce, the quantum 
effects of gravity are no longer negligible, and a description in terms of a classical spacetime 
metric must be abandoned. One needs to learn how quantum helds propagate in a quantum 
gravitational background. 

The quantum theory of cosmological perturbations in the Planck regime was developed 
in [6], based on previous work by Ashtekar, Kaminski and Lewandowski [27], and the result 
can be summarized as follows. In the regime in which perturbations can be treated as test 
helds, i.e., their back-reaction on the FLRW spacetime is small, the equations of motion 
of the operators representing scalar and tensor perturbations are formally the same as the 
equations appearing in classical spacetimes, which in Fourier space are 

Qliv) + 2^{fi) Qh{fi) + {k‘^ +Uiv))Qi:iv) =0, (2.6) 

CL 

77"(r)) + 2^(r))77'(r)) + A;^rg(r)) = 0. (2.7) 

A prime indicates a derivative with respect to conformal time fj, and k = \k\. However, the 
background quantities that appear in these equations, namely a, U and fj, are different from 
their classical counterparts.^ In contrast, a(p), U{ri) and fj are now obtained as complicated 
expectation values in the background quantum state To(a, 0), which involve the associated 
quantum operators d and U as well as the Hamiltonian background operator (see section H.C 
of [7]). Their explicit form will not be needed for this summary. Remarkably, because of the 
formal analogy with quantum held theory in classical spacetimes, it is possible to import all 
the well-developed mathematical machinery from that framework—e.g., Fock quantization, 
transition amplitudes, renormalization and regularization—to build a well-dehned quantum 
held theory in quantum spacetimes. Therefore, once the quantum backgrounds quantities 
a, U have been computed, the evolution of perturbations and computation of observables 
closely follows the formalism commonly used in semiclassical cosmology. All the difficulties 
of the quantum spacetime To(a, 0) are encoded in the expectation values d{fj) and U{fj). 
This formalism goes under the name dressed geometry approach, because the evolution of 
perturbations turns out to be mathematically equivalent to the evolution in a curved space- 
time whose metric has been dressed by quantum gravity ehects. The framework has been 
applied in [35] to explore the phenomenology of quantum states ^ 0 ( 0 , 0 ) which have large 
huctuations in the Planck regime and are therefore not semiclassical in any sense. 

But a signihcant simplihcation to this process appears when \['o(a, 0) is chosen to be 
a highly peaked state whose expectation values are well approximated by the solutions of 


is more convenient to use the gauge-invariant variable Q in place of TZ, and convert its power spectrum 
Pq{x) to the desired power spectrum Pjzix) at the end of inflation, than it is to compute the power 
spectrum of the field TZ{x) directly. 

^ The classical potential that appears in the equation for Qj;{y) \-sU{p) = a^[F(0) r — 210(0)-v/r + 
where r = 3a^0'^ F(0) is the inflaton potential, and F^(0) = dF(0)/d0. 
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the effective equations (2.1) and (2.2). For such highly peaked states, higher moments of 
'l/o(0; 0) are well approximated by powers of the simplest expectation values. Consequently, 
explicit computations show that the scale factor d^fj) appearing in (2.6) and (2.7) reduces 
to the solution a{r]) of the effective equations, and the potential V({r]) takes the same form 
as in the classical theory where the time evolution is now dictated by the effective equations 
of LQC. Therefore, for states Tq that are highly peaked, the computation of power spectra 
and other physical relevant quantities in LQC follows the same steps as in general relativity, 
with the only difference being that the evolution of the scale factor a{r]) and the background 
inflaton held (j){r]) is now replaced by the solution of the effective equations of LQC. 

The quantization of scalar and tensor perturbations now proceeds in the standard way, 
which we summarize here. For scalar perturbations, one begins by decomposing the held 
operator Q(x) as 

/ d^k - - 7 r - -u. 1 - 

^Sf(r,)e'« = J ^ + (2.8) 

where the functions qk{v)^ labeled by k, form a orthogonal basis of the subspace of “positive 
frequency”—more precisely, positive norm—complex solutions to the equation of motion 
(2.6). If the basis elements are chosen to satisfy the normalization condition 

QkivWkiv) - QlivWkiv) = , (2.9) 

a[riY 

then and At satisfy the usual algebra of creation and annihilation operators [Ag, A^,] = 

lit, it] = 0, |ij, 4] = h(27,)H(k- i). 

The vacuum is the state annihilated by all Ag, and the symmetric Fock space is the 
Hilbert space generated by repeatedly operating on this vacuum with creation operators. 
Eq. (2.8) is then the representation of the operator Q in this Hilbert space. It is important 
to emphasize that the dehnition of vacuum is tailored to the choice of “positive-frequency” 
basis functions qk{v)- maximally symmetric backgrounds such as Minkowski or de Sitter, 
one can use the spacetime isometries, together with suitable regularity conditions, to single 
out a preferred basis qk{v) a preferred vacuum state. But in spacetimes with fewer 
isometries, e.g., homogeneous and isotropic FLRW backgrounds with arbitrary scale factor 
a{r]), there is no canonical vacuum, and consequently the notion of particle is ambiguous. 
Note that by using a basis of Fourier modes qkij}) fhat only depend on the length k of the 
wave vector, rather than its direction, one is already restricting to a family of vacua in which 
all members are isotropic and homogeneous; this is manifest in the form of the two-point 
function written below. But there still remains inhnite freedom in the choice of vacuum, 
even within this family.^ 

In a free theory, all vacuum correlation functions can be written in terms of the two-point 
function, which—given a choice of basis qkid )—can be computed as 

/ d^k - 

(lkivi)Q*kiV2) . (2.10) 

^ The adiabatic condition (see [6] for a summary) forces the mode functions qkiv) to approach Minkowski 
positive-frequency modes at a specific rate when the physical momentum k/a{t) is large compared to the 
spacetime scalar curvature. But this is an asymptotic condition, and therefore there are infinitely many 
choices of basis modes qkiq) satisfying it. 
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The relevant observable in cosmology is the two-point fnnction in momentnm space at con¬ 
current times. For a homogeneous and isotropic vacuum, this two-point function is diagonal 
in the two momenta involved, and a simple computation shows it is given by 

^ ^ 27r2 

(0|Qg(r/)Qg,(77)|0) = {2'Kf5{k + k') -j^PQ{k,ri), (2.11) 

where Pq is the power spectrum, which is written in terms of the mode functions as 
Pe(k,ri) = RV, |gfc(? 7 )p. Although this power spectrum of the Mukhanov-Sasaki variable is 
easier to compute, the quantity which is more directly related to observations and thus more 
interesting in inflationary cosmology is the power spectrum of comoving curvature pertur¬ 
bations Pti- Evaluated at the end of inflation r^end, this power spectrum can be obtained 
from Pq as 

PAk) = . (2,12) 

\ <P{Vend) / 

During the inflationary era, sometimes it is useful to write the mode functions qkiv) 
obtained from the pre-inflationary evolution in terms of the modes that dehne the 

Bunch-Davies vacuum during slow-roll inflation 

qkiv) = ^kqk^iv) + Pkqf^*{v), (2.13) 


where qk^{ri) = Hj^\—kri), with (x) a Hankel function of the first kind of 

order /x = ‘i/2 + 2e + 5. Then, Pnik) can be written as the power spectrum for the modes 
times a factor encoding the pre-inflationary evolution which involves the Bogoliubov 
coefficients ak and 13k- 


Pn{k) = P'TlHk)\ai, + h?, 


(2,14) 


where 



/ FT (y^gnd) 
V 0(hend) 



AtiG (H{r,k)V 

f-irik) V 2ir y 


(2,15) 


with the Hubble exit time for the mode k, rjk, dehned by the relation k/a{rjk) = H{rjk). 

The analysis of tensor perturbations is analogous. The held operator is now expanded in 
terms of the elements ek{f]) of a basis of the space of “positive-frequency” complex solutions 
of Eq. (2.7) as 


f{x) = 


d?k 


TM = 


d?k 

(^ 


I3^ek{ri) + IBKeliri) 


Ak-x 


(2.16) 


where the basis functions are normalized to 


efc(h)4*(^) - 4(^)4(^) = 327rG—^ . (2.17) 

a[riY 

The power spectrum for each polarization is given by 

PtA) = h—AAend)? = P?{k) Al + Pl\\ (2.18) 

where P^\k) = h327iG(j^AP^ , and a]/ and I3'l are the Bogoliubov coefficients relating 

CkA) and the Bunch-Davies modes Hi^\—kr]), with u = 3/2 -|- e, during 

slow-roll inhation. 
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C. The LQC power spectrum 

The scalar and tensor power spectra in LQC were computed and analyzed in detail in 
[7] following the theoretical framework we have just described above (see [8-24] for other 
approaches within LQC). Here, we will conclude this section by summarizing some results 
found in [7] and also presenting some new features. 

We have seen that a unique evolution corresponds to each choice of the parameters 
(0B,'nr). For dehniteness, we will consider in this section the power spectra generated by 
choosing m = 1.3 x 10“®, which corresponds to the value that is commonly used in standard 
inflation, and = 1. (Recall that we use Planck units throughout this paper.) We also 
must supply concrete initial conditions for the perturbations. As an illustrative example for 
this section we choose the preferred instantaneous vaeuum introduced in [40]; we will impose 
this vacuum at an initial time 25000 Planck seconds before the bounce, when all the modes 
of interest are inside the curvature radius. As previously discussed, and as will be shown 
explicitly in section IV, other reasonable choices of initial vacua and initial times produce 
power spectra which are all very similar. 

It is useful to provide a qualitative understanding of the physical evolution of pertur¬ 
bations across the bounce. There are two relevant energy scales in the problem. On the 
one hand, LQC introduces a new energy scale /i:LQc/o(tB) = a/Rb/ 6 ~ 3.21 that is directly 
related to the spacetime scalar curvature at the bounce R-q = dSvrpmax ~ 62. (This is the 
maximum value that the curvature attains along its evolution.) A second scale is provided 
by the value of the scalar curvature at the onset of inflation ki/a{ti) = a/Ri/ 6 ~ 10“^. For 
kinetic-dominated bounces, the onset of inflation takes place at ti ~ 5x 10"^ after the bounce.® 
To understand the qualitative features of the pre-inflationary evolution of a Fourier mode 
with co-moving wavenumber k, it is convenient to divide the discussion in three different 
groups: 

• Fourier modes with co-moving wavenumber k > /clqc- These modes are “inside” the 
curvature radius at the bounce (i.e., their wavelengths are smaller than the radius of 
curvature at that time), and will continue to be so until the slow-roll inflationary era. 
Consequently, spacetime curvature will not affect their evolution until the inflationary 
era. One expects that these modes will reach the onset of inflation in the Bunch- 
Davies vacuum, and the hnal power spectrum will have negligible contributions from 
the LQC pre-inflationary evolution. 

• Modes with co-moving wavenumber ki < k < A;lqc- They are “outside” the curvature 
radius at the bounce. But evolution will bring them inside soon after the bounce, 
and they will exit again during slow-roll inflation. This process of “curvature radius¬ 
crossing” enhances the amplitude of the perturbation. In the semi-heuristic language of 
particle creation, the evolution will create quanta as a consequence of the interaction 
with spacetime curvature, and the onset of inflation is reached in an excited state. 
Therefore, we expect the pre-inflationary evolution to affect signihcantly the power 
spectra of those perturbations. 

• Modes with co-moving wavenumber k < ki. They are “outside” the curvature radius 


® Note that this time is different from the onset of slow-roll, ton, dehned previously. Onset of inflation here 
is defined as the beginning of the phase of accelerated expansion, at which e becomes smaller than unity. 
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FIG. 2; The LQC scalar power spectrum for parameter values (j)-Q = 1, m = 1.3 x 10“®, and 
preferred instantaneous vacuum [40] initial data for perturbations at initial time t = —50000. The 
numerically evolved spectrum, shown in gray, is rapidly oscillatory; its average, shown in black, 
has an amplitude which is amplified with respect to the standard predictions of slow-roll inflation 
for modes ki < k < /clqc but agrees with them otherwise. 


at the bounce and will continue to be so all the way until the end of inflation. Since 
these modes do not cross the curvature radius, it is expected that their power spectrum 
remain small as compared to other modes that do cross it. 

Therefore, we expect LQC corrections to be relevant for Fourier modes with co-moving 
wave number ki < k < /clqc- The important question is whether those modes are observ¬ 
able today. To answer this we need to know what physical scales the two scales fengc and 
ki correspond to at the present time. The wavenumbers that we can observe in the CMB 
lie approximately in the interval (/c*/8.9, lOO/c*), where fc*/ao = 0.002 Mpc“^ is the “pivot” 
mode used by the Planck satellite team in parameterizing the primordial power spectrum 
[28], and corresponds to approximately £ = 27 in the angular power spectrum today.^ Be¬ 
cause physical momenta k/a{t) exhibit redshift as the universe expands, the values of the 
quantities ^Lgc/ao and ki/ao depend on how much expansion has occurred from the bounce 
and the onset of inflation, respectively, until the present time. The amount of expansion, in 
turn, is controlled in part by the values chosen for 0 b and m. Thus the question of whether 
or not the scale /clqc is observable today, and therefore whether the effects of LQC physics 
may be potentially imprinted in the CMB, will depend on the value of the parameters 0 b 
and m. 

The above discussion, although very useful to create an intuitive picture, is heuristic. 
It needs to be supported by explicit computations. Figures 2 and 3 show the scalar and 
tensor power spectrum, respectively, obtained after numerically evolving a large range of 


^ -^ki^/ao equals the Hubble rate today. 
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FIG. 3: The LQC tensor power spectrum for parameter values (/>b = 1, m = 1.3 x 10“®, and 
preferred instantaneous vacuum [40] initial data for perturbations at initial time t = —50000. The 
numerically evolved spectrum is shown in gray, and its average in black. The qualitative behaviour 
is similar to the scalar spectrum shown in Fig. 2. 


co-moving Fourier modes fc.® As previously mentioned, this figure corresponds to the back¬ 
ground values 0B = 1, = 1-3 x 10“®, and preferred instantaneous vacuum initial data [40] 

for perturbations specihed in the pre-bounce, contracting phase at f = —50000, when all 
interesting modes are well inside the curvature radius. The figure clearly shows the three 
regions in fc-space previously specihed. For k > /clqc the power spectrum has negligible 
LQC corrections, and therefore the results agree with the inhationary prediction obtained 
using Bunch-Davies vacuum initial data at the onset of inhation. In contrast, the inter¬ 
mediate region ki < k < /clqc shows a signihcant enhancement coming from the LQC 
pre-inhationary evolution. Finally, the power spectrum for the longest wavelengths k < ki 
is largely suppressed. 

For the chosen values 0 b = 1 and m = 1.3 x 10“®, the LQC scale today /cLQc/no is around 
one-third of the pivot scale fc*/ao = 0.002 Mpc“^. Therefore, the CMB would show some 
LQC contributions, although only on the largest angular scales E < 15. The scale ki/ao 
is approximately 10“^ times k^/ao] equivalently, it corresponds to a wavelength 10^ times 
larger than the Hubble radius today. 

To summarize, the LQC corrections to the primordial spectrum of cosmic perturbations 
are more pronounced for low values of k (i.e., long wavelengths). Depending on the values of 


® The behavior of the numerically evolved scalar power spectrum is highly oscillatory in k\ this phenomenon 
is generically expected for departures from the standard inflationary paradigm (see e.g. [41-43]). All state¬ 
ments we make about the power spectrum should be interpreted as referring to its value after averaging 
over these high-frequency oscillations, which are so rapid in k as to be unmeasurable by any realistic 
observation. 
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the background parameters 0 b and m, it is possible that LQC corrections begin to appear 
for the longest-wavelength modes that we can observe, or only for wavelengths that are 
larger than the Hubble radius today and thus not directly observable. The power spectrum 
is further amplihed for yet larger super-Hubble wavelengths. But remarkably, the power 
does not grow unboundedly for smaller k; on the contrary, it reaches a maximum around 
k = kj and then decreases quite abruptly for lower k values. 

The power enhancement of super-Hubble modes in the range between k\ and /clqc is 
quite interesting and, although counter-intuitive at first, it may lead to additional observable 
effects. This could be the case if strong correlations between observable and super-Hubble 
modes happen to exist, i.e., if there is considerable non-Gaussianity between these two sets 
of modes. Under this circumstances, super-Hubble scales can indeed affect the observed 
power spectrum, introducing a modulation superimposed on the power spectrum we have 
showed in Fig. 2. Such modulation could help us to understand the origin of some anomalies 
discovered by WMAP and conhrmed by Planck at the largest angular scales of the CMB. The 
computations of the non-Gaussianity arising as a consequence of the large power spectrum in 
super-Hubble modes, and its effects on large angular scales in the GMB, has been analyzed 
in [44]. 


III. EXPLORING THE PARAMETER SPACE 

The previous section summarized the theoretical framework needed to propagate cosmic 
perturbations in the early universe and to compute the power spectra in LQG for a specihed 
choice of initial parameters 0 b and m. In this section, we consider pairs (0b, m) through¬ 
out the relevant parameter space dehned in Eq. (2.4), and inquire for which such pairs the 
resulting power spectra are consistent with the strict constraints coming from observations. 
Then, among this region in the parameter space which is observationally viable, we addi¬ 
tionally identify the subset of those points which nonetheless incorporate signihcant LQG 
modifications, and we characterize what those modihcations are and how they might be 
distinguished in future observations. 


A. Observational constraints 

Our ultimate goal is to contrast the result of our computations with observations. To that 
end, in this section we summarize the constraints stemming from GMB observations that are 
relevant for our computations. We will use Planck 2013 results [28] for most of our analysis. 
Observational constraints come entirely from scalar perturbations; for tensor perturbations 
we have only an upper bound on their amplitude, coming from a joint Planck-BICEP2/Keck 
Array analysis [45]: r(A;*) < 0.12 (95% GL). 

1. Amplitude and spectral index of the scalar power spectrum. By using a phenomeno¬ 
logical parametrization of the primordial power spectrum given by P-jiik) = Ag 
the Planck data, in combination with WMAP and BAO, provide the following values for 
the best ht of the scalar amplitude Ag and spectral index Ug [28]: 

A, = (2.196+°;°!) X 10"® > = 0.9643 ± 0.0059 . (3.1) 

where la uncertainty ranges are indicated. 
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2. Running of the spectral index. When a running as = dns/dln/c is included in the 
parametrization, Pn{k) = Ag the Planck 2013 data produce 


a, = -0.0013 ±0.009, 


(3.2) 


again displaying a la uncertainty range. This running is compatible with zero and even a 
positive value at a 1.5a level. The inclusion of running does not improve significantly the 
maximum likelihood of the parametrization. 

The recently released new 2015 Planck data [29] provide slightly different values for Ag, 
Ug and ag. The impact this change produces on our conclusions is negligible. 


3. Number of e-folds of inflation. The number of e-folds A/”* between the time tk^ at 
which the pivot scale /c* left the Hubble radius during inflation and the end of the inflationary 
era fend is constrained [46]. (Recall that fend is dehned here as the time when the slow-roll 
parameter e = reaches the value e = 1 for the hrst time after slow-roll.) The origin 

of this constraint is as follows: The pivot scale /c* is, at the present time, 8.9 times smaller 
than the Hubble scale, i.e., Hq = 8.9 k^^/aQ. On the other hand, the Hubble exit time of the 
mode /c* during inflation, tk^, is dehned by the relation H(tkfl = k.^/a(tkfl. Therefore, by 
eliminating /c* from these two equations we have a relation between quantities at f^^ and at 
the present time: iJoOo/8-9 = H(tkfl a(tkfl. From this expression it is straightforward to 
obtain 


A/”* = In 


® (fend) 

a(ffcj 


In 8.9 ± In 


/ Ci(fend) Ff(ffc^) 

V Hq 


(3.3) 


If we make the extreme assumption that the process of reheating is instantaneous and non- 
dissipative (meaning all energy in the inhaton potential is converted into radiation), then, 
for the quadratic potential that we use in this paper. A/”* can be estimated as approximately 
61. A more realistic reheating process could both increase and decrease this quantity; taking 
into account this uncertainty, we allow the conservative range 


50 < A/; < 70 . 


(3.4) 


B. Constraining the parameter space 

In this section we report the main computational results of this paper. We have written 
a numerical code which systematically computes the scalar and tensor power spectra for 
values of and m in the relevant parameter space discussed in section H—pairs {(f)B,rn) 
with 0 < ± 0.90—and contrasts the result against observations. The code is writ¬ 

ten in Mathematical it evolves the mode functions using Mathematical native numerical 
differential-equation solver, which uses an LSODA approach, switching between a non-stiff 
Adams method and a stiff Gear backward differentiation formula method. The code was run 
on the “Philip” high-performance computer cluster at Louisiana State University, where it 
requires approximately 10 processor-hours per point in the parameter space. 

Using the observational constraints spelled out in the previous subsection, the code clas- 
sihes pairs (^b,’^) into the following three (overlapping) categories: (i) Pairs for which 
the scalar power spectrum is compatible with Planck observations for Ag and Ug. (2) The 
subset of the those points for which the scalar or tensor power spectra contain signihcant 
LQC contributions. (3) Points for which the number of e-folds Afl satisfies the constraint 
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spelled out at the end of previous subsection: 50 < A/”* < 70. In the following we provide 
further explanation of each category and show the results. 

For dehniteness, the computations in this section are done using the preferred instanta¬ 
neous vacuum initial condition for perturbations imposed at t = —0.2.^ Section IV discusses 
the use of other initial conditions and shows that the results obtained here are unchanged 
for different reasonable choices of initial vacuum or initial time. 


(1) Points (0B, m) compatible with Planck observations for Ag and Ug. Pairs (0 b, m) 
belong to this category if the scalar power spectrum contains at least one value of k for which 
the amplitude and tilt agree with Eq. (3.1) inside their joint la uncertainty region.(The 
uncertainty ranges make it possible for more than one k to satisfy this condition.) These k 
values are candidates for the pivot scale A;*. 

Points in the plane (0 b, ''^) satisfying this condition are in the region of Fig. 4 outlined 
by the thick black line, which corresponds to 0 b ^ 0.8 and 0.9 x 10“® ^ ^ 1-6 x 10“®. 

Note that these points cover a small portion of the theoretically allowed parameter space. 
However, that is already the case in the standard inflationary paradigm without LQC, or 
more generaly in any physical model after contrasting with observational data. 

It may seem surprising at hrst that it suffices to check the values of Ag and Ug at the 
single mode /c* to ensure compatibility with observations. The reason is that, because the 
LQC corrections are more important for low /c’s, it is guaranteed that pairs (0b, m) with 
appropriate amplitude and tilt at /c* also show negligible deviation from a power law for all 
k > /c*. Therefore they are compatible with Planck observations. 


(2) Points (0B, m) compatible with Planck observations for Ag and Ug that incorporate 
significant corrections from LQC. This category is made of the subset of points from cat¬ 
egory (i) which show at least a 10% contribution from LQC-physics. To more precisely 
dehne what “contribution from LQC-physics” means, it is convenient to write the scalar 
power spectrum in terms of the auxiliary power spectrum dehned in Eq. (2.15); then, the 
amplitude and tilt can be written as 


Pn{k)=p!^\k)\ak + Pk\\ 


ng{k) — 1 = nf'\k) — 1 -|- 


dln|Q;fc + 0fcP 

d In /c 


(3.5) 


where p!^\k) was given in Eq. (2.15) and nf\k) — 1 = —Ae{rik) — 25{rjk). These expressions 
are very useful because they neatly codify the contribution of pre-inflationary physics into the 
Bogoliubov coefficients au and 0^, while the standard inflationary contributions are included 
in P^\k) and nf\k). Therefore, we will say that a given couple (0b, rn) contains signihcant 
contributions from LQC physics when the ak and (dk contributions modify somewhere (in 
/c-space) the inflationary results by at least 10%. The factor \ak + 0a:P can be rewritten, by 


® We use t = —0.2 rather than the bounce t = 0 because the preferred instantaneous vacuum initial data 
are not well defined at the bounce time. 

Note that the Icr region in the {As,ns) plane is not simply the Cartesian product of the individual Icr 
regions given in Eq. (3.1) for Ag and Ug. Rather, we model Planck's measured values for the magnitudes 
and uncertainties of Ag and Ug as following skew-normal and normal, respectively, distributions, and find 
the (roughly elliptical) contour in the (H^, ng) plane for which the joint probability is no greater than la. 
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FIG. 4; Regions of the parameter space that yield LQC power spectra meeting various 

criteria. The thick black outline demarcates the points (i) consistent with the Planck 2013 ob¬ 
servational constraints on the amplitude and tilt of the scalar power spectrum. The gray region 
{2) indicates the subset of such points for which the spectra contain potentially observable LQC 
contributions. Finally, the striped region (5) indicates the points consistent with the constraint 
on the number of e-folds A/"*. Several values at selected points from this hgure (indicated by the 
black shapes) are given in Table I. Due to prohibitive constraints on computational resources, the 
shapes of the regions (i) and (5) have been extrapolated for the largest shown values of (pB- 



m 

4>b 


r{k^) 

nt{k^) 

r/nt 

as{k^) 



T 

0.63 X 10"® 

2.12 

8.3 

0.07 

-0.035 

-2.0 

2.5 X 10-3 

5.49 X 10-® 

0.0044 

• 

0.96 X 10"® 

1.37 

10.7 

0.11 

-0.018 

-6.0 

3.3 X 10-4 

6.70 X 10-® 

0.0068 

■ 

1.06 X 10"® 

1.22 

10.4 

0.12 

-0.021 

-5.8 

4.3 X 10-4 

7.03 X 10-® 

0.0075 

♦ 

1.28 X 10“® 

1.00 

16.0 

0.14 

-0.018 

-7.9 

2.9 X 10-® 

7.81 X 10-® 

0.0089 

♦ 

1.51 X 10“® 

0.79 

15.0 

0.17 

-0.021 

-7.9 

8.2 X 10-® 

8.46 X 10-® 

0.0106 


TABLE 1; Values of the tensor-to-scalar ratio r, the tensor spectral index nt, the quotient r/nt 
appearing in the consistency relation, the running of the scalar spectral index ag, and the value of 
the Hubble rate H and the slow-roll parameter e at rj^* for various points selected from 

the parameter space. The symbols appearing at the left of the table correspond to the symbols in 
Fig. 4 above. 
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taking into account the normalization condition la/cp — l/dfcP = 1, as H-2|/3fcp + 2 Re 
We observe in our computations that the interference term Reja^/d^} is, for the range of 
/c’s relevant for this section, highly oscillatory with zero average. Therefore, we average it 
out and will say that a given couple contains signihcant contributions from LQC 

physics when 


2\f3k\^ > 0.1 or 


1 dln|l + 2/jfcp 
nf\k)-l dln/c 


(3.6) 


for some value of k in the observable range. 

The pairs (0 b, m) satisfying this condition appear in gray in Fig. 4. They are distributed 
narrowly around a curve given approximately by 0 b = 1-3 x 10“®/m. For all these points, 
the LQC corrections appear for the lowest values of k. 

Table I indicates the values of the tensor-to-scalar ratio r, tensor spectral index n*, scalar 
running o;*, and other quantities of interest for some points in this category. We observe 
that: 


• The LQC corrections increase for lower values of m along the curve 0 b = 1.3 x 10“®/m. 
In previous analysis [7] it was pointed out that the LQC corrections decrease for larger 
values of 0 b when m is held hxed. This is manifest in Fig. 4 where, for any given value 
of m in the range (0.96,1.6) x 10“®, LQC corrections weaken when we move vertically 
upward in the hgure (becoming negligible once out of the gray region). When the 
freedom in m is included, our results indicate that LQC corrections increase rapidly 
for lower values of m, more than compensating for the effect of increasing 0 b, when 
we approximately follow the curve 0 b = 1.3 x 10“®/m toward smaller m. 

• The trend of the LQC-corrections is to make n* more negative, and to increase 
This can be understood by simple inspection of Figs. 2 and 3, where the corrections 
increase for lower k, therefore making the spectral index more negative and increasing 
its running. 

• LQC corrections decrease both the value of the slow-roll parameter and the Hubble 
rate at the time the reference mode /c* exited the Hubble radius during inflation. The 
reason these values are modihed is as follows. The scalar amplitude Ag and spectral 
index Ug —whose values are fixed, up to error bars, by the observational constraints of 
Eq. (3.1)—are given in LQC by^^ 


A. 


Pn{k,) = h 


AnG 



2 

+ 0A:* P 


(3.7) 


and 


ng{k^) = 1 - 4e(?7fcJ + 


din \ak + 


dln/c 


fc=fc* 


(3.8) 


In standard inflation with vacuum initial conditions, the last term in Eq. (3.8) vanishes, 
while it is negative in LQC; to maintain the observationally mandated value of Ug, the 


See Eqs. (3.5). The latter of these equations has been particularized here for the quadratic potential 
discussed in this paper; however, the same argument holds for other forms of V (d) ■ 
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value of e{r]k^) must be smaller than in standard inflation. Then, since Ittfc + /dfcP > 1, 
Eq. (3.7) implies that LQC corrections make H{r]kJ smaller as well. 

As a consequence of the decreased value of Hirjk^) in LQC, the energy scale of inflation 
at the time observable perturbations were generated is reduced. 

• LQC corrections tend to decrease r. As pointed out in [7], the tensor-to-scalar ratio 
in LQC is 

= ^ = (3,9) 

where we have made use of Eqs. (2.14) and (2.18), and in the last equality we have used 
the fact that the Bogoliubov coefficients for scalar and tensor perturbations are very 
similar. This expression looks exactly the same as the result that one would obtain 
in standard inflation with Bunch-Davies vacuum initial conditions (i.e., without LQC 
corrections). But we noted above that the corrections decrease the value of e(? 7 fc^); 
thus they in turn decrease the predicted value of r(A;*). 

Therefore, LQC helps to alleviate the observational constraint on the potential. 
However, when all observational constraints are imposed, the corrections on r are 
small, and this potential still remains close to the border of the 95% CL region. We 
will have to wait to polarization to know if the potential is dehnitively ruled 

out. But we emphasize that even if that turns out to be the case, the results of this 
paper will remain valid for other potentials (e.g., see [47] for the detailed analysis of 
the Starobinsky potential in LQC). 

• The ratio r/rit can depart from the standard consistency relation of inflation r/rit = 
—8; the tendency is for this ratio to become less negative with decreasing m (see Table 
I). More precisely, in LQC the consistency relation becomes [7] 


r = —8nt + 


dln/c 


(3.10) 


The extra term ^ ^ is negative in LQC, therefore adding a positive term to 

the ratio r/ut (recall Ut is negative). 


(3) Points (0B, m) compatible with constraint on A/”*. In this category we add the con¬ 
ditions on the number of e-folds to the points in category 1. As described previously in 
subsection III A, consistency with the present size of /i;*/ao requires k^, to have exited the 
Hubble radius between 50 and 70 e-folds before the end of inflation, 50 < A/”* <70. As ex¬ 
plained above, imposing that the amplitude and tilt fall within the joint la region observed 
by Planck provides a range of candidate values of A;* at each point in category 1. We have 
marked with stripes in Fig. 4 the region of points for which at least one of these candidates 
for /c* also satishes the constraint on A/”*. 


IV. INITIAL CONDITIONS FOR PERTURBATIONS 

To compute the power spectrum, in addition to specifying the background parameters 
0B and m, one also needs to provide initial conditions for perturbations. This is done 
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by specifying the quantum state of perturbations at some instant of time, at least for the 
Fourier modes of observational interest. At the practical level, as explained in section II, 
this is achieved by providing initial data for the modes functions qk{v) and 6 ^( 77 ) of scalar 
and tensor modes, respectively, for the values of k we are interested in. As also discussed 
in that section, in a generic FLRW spacetime there is no preferred or canonical choice for 
these initial data. The first choice one needs to make is when to specify initial data. Two 
natural times to impose “vacuum” initial conditions are the bounce and the far past (see 
section III C of [7] for physical arguments in favor of imposing vacuum initial conditions 
at the bounce in LQC). But even after a choice of initial time is made, there is still the 
ambiguity of which “vacuum” we choose. In this section, we analyze the sensitivity of the 
observable predictions described in section III to different choices. The conclusion will be 
that different reasonable choices of adiabatic vacua, imposed at the bounce or at some time 
prior, all produce very similar results for observable modes. Thus the main conclusions of 
section III are robust. 


A. Sensitivity to the choice of vacuum initial conditions at a given time 

We have explored the sensitivity of observable quantities (power spectra and spectral in¬ 
dices) to different choices of vacuum state. At the time of the bounce, our analysis reduces 
to the one in [7], and our conclusions are the same: observable quantities are quite insen¬ 
sitive to the particular choice of vacuum at the bounce time. We therefore refer the reader 
to [7] for further details. It is important—although perhaps obvious—to emphasize that 
these conclusions cannot be extrapolated to any vacuum state one can write. An arbitrary 
Bogoliubov transformation of a given adiabatic vacuum, with appropriate fall-off conditions 
as /c —)■ cx), is also a legitimate vacuum. Then, by choosing Bogoliubov coefficients 13k ap¬ 
propriately, one can modify the observable quantities arbitrarily much. Our conclusions, as 
those in [7], are therefore restricted to the small set of states we have explored, which have 
been selected by demanding extra physical conditions that make those states reasonable 
candidates for the ground state (see e.g. [40]). 

We have extended the analysis in [7] to times before the bounce, obtaining similar con¬ 
clusions. As an example. Fig. 5 shows the scalar power spectrum for two states dehned 
at 50000 Planck seconds before the bounce, namely the preferred instantaneous vacuum 
introduced in [40] and the state with Minkowski-vacuum-like^^ initial data at that time. 
At times sufficiently before the bounce, all observable modes had physical momentum well 
above the spacetime curvature scale (i.e., the observable modes were well “inside” the cur¬ 
vature radius). In that limit, all vacua of at least second adiabatic order differ only by terms 
of order where R{ri) is the Ricci scalar. Therefore, in that regime different choices 

of adiabatic vacua all produce very similar result. In the extreme limit 77 —)■ — 00 there is a 
preferred notion of ground state: the Minokoski vacuum—or the Bunch-Davies vacuum if a 
positive cosmological constant is included in the model. 


This state has initial data at y = rjo given by qk(j]o) = ^/[a(r]o)V^] and q'kijlo) = —ikqkiyo) for all 
k. Strictly speaking the resulting state is not adiabatic (it is only zeroth-order adiabatic). However, one 
can modify these initial data for values of k much larger that the largest k we can observe in order to 
make the resulting state fourth-order adiabatic; observable quantities are unaffected by this process of 
“adiabatization”. 
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FIG. 5; Averaged LQC scalar power spectrum for two different choices of initial data for pertur¬ 
bations specified at 50000 Planck seconds before the bounce. The two spectra are nearly indistin¬ 
guishable for observable modes. The background parameters used for both spectra in the figure 
are m = 1.3 x 10“® and cpB = 1- 



FIG. 6: Averaged LQG scalar power spectrum for Minkowski-like initial data specified at different 
times, t = 0 corresponds to the bounce, and values of time are given in Planck units. The power 
spectrum changes significantly for very low values of k, indicating that those modes are sensitive 
to the pre-bounce evolution. On the other hand, all curves are very similar for observable modes, 
except for the lowest observable values of k. The background parameters used for both spectra in 
the figure are m = 1.3 x 10“® and cpB = 1- 
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B. Sensitivity to the choice of time at which vacuum initial conditions are imposed 

We have also explored the extent to which observable quantities are affected by imposing 
the same notion of vacuum initial data for perturbations at different times. Again, we 
hud that predictions are quite insensitive to the specific time chosen for observable modes. 
However, we also observe that the power spectrum for modes with smaller wavenumber k 
(longer wavelength), which today are super-Hubble, are significantly affected. Our result are 
illustrated in Fig. 6, which shows the scalar power spectra for 0 b = 1 and m = 1.3 x 10“® 
arising from Minkowski-like initial conditions at different times. The hgure shows that the 
power spectrum for modes k with small ratio k/k* changes signihcantly when initial data 
is specihed at different times. However, the power spectrum is unchanged for modes with 
k/ki, > 1. The physical reason is that modes with low k exit the curvature radius well 
before the bounce and consequently their pre-bounce evolution is significantly affected by 
the spacetime curvature. On the contrary, modes with /c/Zc* ~ 1 only “feel” the spacetime 
curvature very close to the bounce time, and therefore their power spectrum is insensitive 
to which time we choose to specify initial data in the contracting phase. 


V. CONCLUSIONS 

In this work we have used observational data to constrain the parameter space of the phe¬ 
nomenological sector of loop quantum cosmology, and we have provided a detailed analysis 
of the shape of quantum gravitational corrections to observable quantities along this param¬ 
eter space. We have emphasized that the freedoms in the value of the parameters appearing 
in the inflaton potential V (0) are not hxed by observations alone, and must be included in 
the parameter space. The computations presented here require the use of high-performance 
computing. It is our view that this is a necessary task to have full control of the phenomeno¬ 
logically allowed range of parameters and the predictions of the model thereon. We expect 
our results to be particularly useful for contrasting the theory with the forthcoming data on 
CMB polarization, as well as to extend the phenomenological explorations in LQC beyond 
the power spectrum—as for instance in the computations of non-Gaussianity [44]. 

We hnd particularly remarkable the tight constraints that current observations impose 
on the parameter space of quantum gravity. When thinking of quantum gravity corrections 
to observable quantities, one tends to imagine minuscule contributions, hence leaving large 
freedom for speculation without conflicting with observations. Our analysis shows that this 
is far from the case. On the contrary, we have seen that observations strongly constrain 
the potential effects of quantum gravity. This is more clearly illustrated in Fig. 4, where 
LQC corrections are constrained to a tiny subset of points, distributed in an almost one¬ 
dimensional strip, across the theoretically allowed parameter space. 
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